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Abstract 



We studied the leading supersymmetric contribution to top-antitop threshold 
production using the NRQCD framework. The one-loop matching to the po- 
tential and the Wilson coefficient of the leading 3 Si production current were 
considered. We point out that the leading correction to the potential is zero 
due to SU(3) C gauge invariance. This is true in general for any new physics 
that enters above the electroweak scale. The shape of the top quark pair pro- 
duction cross section is therefore almost unaffected near threshold, allowing 
a precise determination of the top quark mass based on the Standard Model 
calculations. The supersymmetric correction to the Wilson coefficient c\ of 
the production current decouples for heavy super particles. Its contribution 
is smaller than the Standard Model next-to-next-leading-log results. 
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It is likely that top quark pair production near threshold will be investigated intensively 
in the future at an e + e~ linear collider. It can be used to determine the top quark mass, 
decay width, top Yukawa coupling and strong coupling constant. The Standard Model (SM) 
QCD prediction has been studied in great detail by several groups |l|J^] in the framework of 
NRQCD (nonrelativistic QCD) [|3] up to next-to-next-leading-log (NNLL) order. Near the 
threshold, the velocity of the top quark v ~ a s ~ 0.15 is small. However, the relevant energy 
scales rrit ~ 175 GeV, m t v ~ 25 GeV, mtv 2 ~ 4GeV, together with top decay width T t ~ 1.5 
GeV are much larger than Aqcd- Therefore, top quark production can be calculated to high 
precision using perturbative methods. 

Schematically, the tt cross section a t i can be expressed as a series: 
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The velocity renormalization group (VRG) |4| summation of a s hav and Coulomb summation 
of (Xg/v are performed to all orders. SM QCD calculations have been performed at NNLL 
order ||, and the remaining theoretical uncertainties are likely to be at the few percent 
level. 

Supersymmetry (SUSY) is the most prominent extension of the SM, and gives new 
contributions to top quark production through loop corrections. Supersymmetric one-loop 
corrections to the process e + e~ — > tt above the threshold region have been studied in ||. 
For scales below m t , heavy super particles decouple and they only affect tt pair produc- 
tion in the threshold region through matching conditions at the scale \i — mt- For heavy 
super particles, the supersymmetric contribution is of the order of a t i ~ va s (m t / 'Msusy)" ', 
where Msusy is the scale of the super particle masses, which is smaller than 1 TeV if su- 
persymmetry plays a role in explaining the hierarchy puzzle. Therefore, the leading SUSY 
correction could be comparable to the SM NNLL corrections. In this paper, we will study 
the supersymmetric contribution to top pair threshold production through one-loop match- 
ing conditions, neglecting subleading corrections which are suppressed by higher powers of 
a s or v. We do not consider SUSY loop corrections to the width of the top. They were 
considered in ||. 

For tt production in e + e~ annihilation via a photon or Z boson exchange, the normalized 
cross section in the full QCD theory (at cm. energy y/s) is related to 



Im 



d 4 xe iqx {<d | TJ{x) J f (0)|0> 



(2) 



where q = (\/s, 0) and J is the production current. In NRQCD, The non-relativistic 3 Si 
current for top quark production near threshold can be written as 0] 

J = Cl O P1 + . . . , (3) 

where O p l = ippCr(ia2)x*—p an d the ellipses denote terms suppressed by powers of v. 
The correlator of the currents J in Eq. (^j is replaced by the correlator of P1 , which 
is proportional to the coordinate space Green function G(0,0), with its momentum space 
Fourier transformation G(p,p') being a solution of the Schrodinger equation ||: 



2 



p 2 

— -E~iT t 
m 



G(p,p') + J D n q^V(p,q)G(q,p') = (2n) n 5^(p - p% (4) 



where m is the heavy quark pole mass, E = (y/s — 2m), n = 3 — 2e and D n q = 
e nE (47T)~ e d n q/(2ir) n . V(p,q) is the color singlet quark-antiquark potential in NRQCD, 
with the leading term being just the Coulomb potential. Therefore, the cross section for 
ti production near threshold is proportional to the product of the coefficient c\ with the 
associated Green function. 

The main ingredients in the calculation of the top pair production cross section using 
NRQCD are the potential V(p, q) for a quark-antiquark pair and the Wilson coefficient 
Ci of the production current, which are obtained via matching to the coefficients of the 
relevant operators calculated at the scale /i = m t , then running down to lower scales using 
the VRG y|. Recent SM calculations at NNLL order j7|,||] are likely to have reduced the 
theoretical uncertainties to a few percent level. In what follows, we will consider the SUSY 
contribution to both the potential for a quark-antiquark pair and the Wilson coefficient c\ 
of the 3 Si current, through one- loop matching conditions. 

The potential interaction in NRQCD is written as 

C = -V a pxr(p,p')^ P ' a ipp f3 X^-p'xX-P T , (5) 

where a, j3, A, r are the color and spin indices. ip p and Xp are m the SU(3) C representation 
of 3 and 3, which annihilate a quark and an antiquark, respectively. The leading term in V 
is I: 

y(T) 

V = {T A ®f A )^ r) (6) 
k 

where k = p' — p. At leading order, the SUSY contribution to the Coulomb term Vc is 
always zero by SU(3) C gauge invariance. Integrating out heavy supersymmetric particles gen- 
erates gauge invariant local operators at the top quark scale. In the transition to NRQCD at 
leading order in a s , only the kinetic term for the top quark generates the Coulomb potential. 
The subleading terms in the potential could get contributions from supersymmetric loops. 
However, they are suppressed by powers of v relative to the leading Coulomb potential. 

For more general types of new physics that enter above the electroweak scale, and con- 
tribute to top quark pair production via perturbative matching at \i = m t , the Coulomb 
potential is also protected due to the requirement of SU(3) C gauge invariance. Therefore, the 
shape of the top pair production cross section, which is related only to the Green function 
in Eq. (Pi), will stay almost unchanged, allowing a precise determination of the top quark 
mass based solely on the SM calculations, with uncertainties that may be as small as 100 
MeV 0. 

The Feynman diagrams that contribute to the leading Wilson coefficient c\ of the 3 S*i 
current are shown in Fig. |I|. The SUSY corrections <5ci 7 or bc\z to the production current 
associated with a photon or Z boson are defined as: 

photon coupling : A(|e)O p ^ciy + 5ci 7 ); (7) 



Z coupling: ZJgf + g\{ \ - § sin 2 9 w )O p x {c xz + 6c lz ), (8) 
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FIG. 1. One-loop supersymmetric correction of the tt pair production current. The arrow 
shows the flow of quark number and q can be either i\ or ti- 



where at leading order in perturbative QCD: ci 7j z = 1. We consider only the contribu- 
tions from stop t and gluino, assuming that the flavor mixing between generations in the 
gluino-quark-squark couplings is small. The contributions from SU(2)x,, U(l)y gauginos are 
suppressed by the electroweak couplings. Let m^ denote the masses of the two physical 
stop mass eigenstates, with the mixing angle 9: 
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and m g denote the gluino mass. The loop corrections to the production current associated 
with a photon or Z boson give (neglecting terms suppressed by higher powers of v) : 



5ci 7 = C F —[F(m h ,m~ g ,m t ) + F(m~ t2 ,m~ g ,mt) + G(m~ tl ,m~ g ,m t ) + G(m i2 ,m g ,m t )\; 
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where 
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and C F = 4/3 for SU(3) C . 

In realistic SUSY models, for the squark mass matrix in the (q^, qn) basis, the diagonal 
matrix elements M 2 LRR ~ ^Iusy; while the off-diagonal matrix element M 2 R ~ MsusY m t- 
We chose the mass parameters randomly in the intervals: 
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(200 GeV) 2 < Ml L RR < (1000 GeV) 2 , 
-(1000GeV)m t < M\ R < (1000GeV)m 4 , 

200GeV< \m~ g \ < 1000GeV, (14) 

with the requirement that the physical stop masses stay within the range of 200 GeV < 
m~ tl < mj 2 < 1000 GeV. The numerical results for the SUSY correction to the Wilson 
coefficient c\ is shown in Fig. || for the currents associated with a photon [(a), (b) and (c)] 
and a Z boson [(d), (e) and (f)], respectively. The corrections relative to the tree level value 
are small in both cases: —0.2% < <5ci 7 < 0.3% with 97% of the points falling between ±0.1% 
for the photon, and —0.6% < 8c\z < 0.6% with 86% of the points falling between ±0.1% for 
the Z. The enhancement in the correction to the Z current is due to the mixing between 
the ti and tn, which alleviates the partial cancellation between the tree level left and right 
handed Z-current. 

5c\ goes to zero for large squark and gluino masses, since the SUSY contributions de- 
couple for heavy super particles. However, the gluino decouples slower, with 5c\ ~ m t /\mg\ 
[unlike 5c\ ~ (m t /mq) 2 for the squark], as can be seen in G(m£., rrig, m t ) in Eq. ([13]). The 
largest contribution comes from the term proportional to sin 26. For the Z current, the con- 
tribution from the other terms are enhanced due to the left-right mixing in the stop sector, 
which causes the shift in the 9 dependence in Fig. 0(f). 

The SUSY corrections are small compared with the 0{a 2 s ) NNLL SM contribution to C\ 
They are also smaller than the theoretical uncertainties in SM predictions at NNLL level 
0, which are about 2 — 3% due to the variation in z/, for soft and ultrasoft renormalization 
scales [is = m t^ and fius = rn t v 2 . Therefore, the corrections caused by the existence of 
heavy supersymmetric particles are negligible and it is sufficient to extract the top quark 
mass, decay width and top Yukawa coupling from the ti pair production near threshold 
based on SM calculations. 

In summary, we calculated the leading SUSY contribution to top-antitop threshold pro- 
duction via its corrections to the matching condition for the potential and the Wilson co- 
efficient Ci of the 3 Si production current associated with a photon or Z boson using the 
NRQCD framework. The correction to the potential is zero (at leading order in v ) due to 
SU(3) C gauge invariance, which is true in general for any new physics that enters above the 
electroweak scale. The shape of the top pair production cross section near threshold is unaf- 
fected, which can therefore lead to a precise determination of the top quark mass based on 
the SM calculations. The correction to the coefficient C\ of the production current is within 
±0.1% in most cases, which is smaller than the SM NNLL corrections and the theoretical 
uncertainties from higher orders in QCD perturbation theory. 

We have only considered the leading SUSY contribution from the gluino and stop, ne- 
glecting terms proportional to higher powers of a s and v. Flavor mixing between gener- 
ations is usually small and not taken into consideration here. Neutralino, chargino and 
other squarks would contribute through electroweak loops, however their corrections are 
suppressed by the weak coupling constant instead of strong coupling. If there are super par- 
ticles lighter than top quark, additional contributions to the running of the quark- ant iquark 
potential and the Wilson coefficient of the production current from mt to lower scales should 
be taken into account. 
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FIG. 2. Scatter plot of SUSY correction to c±, for the currents associated with a photon [Fig. (a), 
(b), (c)] and a Z boson [Fig. (d), (e), (f)]. and nig are the masses of the lighter stop and gluino, 
respectively, and 9 is the stop mixing angle. The range of the SUSY mass parameters are chosen 
as described in the paper. 
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